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The impurity concentration dependence of superconducting transition 

temperature T c is studied within the framework of two band model of order 
parameter anisotropy. Both intraband and interband electron scattering by 
nonmagnetic and magnetic impurities are taken into account. It is demon- 
strated that at various values of model parameters both types of impurity 
concentration dependence of T c are possible: weak dependence typical of the 
model with s-wave order parameter and strong suppression of T c by impuri- 
ties with critical impurity concentration at which T c = 0. It is found that in 
some cases there is a possibility of a rise of critical temperature with increas- 
ing impurity concentration. The obtained results are consistent with existing 
experimental data for the high-temperature superconductors. 
The most important and still unresolved problem in the theory of high-temperature 
superconductivity (HTSC) is the mechanism of formation of electron-electron attraction re- 
sponsible for unusual properties of these compounds. In this connection an investigation of 
order parameter symmetry is of great importance. On the one hand, a variety of experimen- 
tal data suggests strong anisotropy of order parameter |l|J^]. On the other hand, a number 
of experiments displays the behavior typical of models with s-wave order parameter ||f|]. It 
was shown earlier that order parameter anisotropy can be explained by symmetry properties 
of superconductor crystal lattice within the framework of the universal multi-band model 
0. The approach offers a non-contradictory explanation of a number of experimental data 
for high-temperature superconductors irrespective of the nature of superconducting pairing. 

Order parameter symmetry proves to be essential if the effect of impurities on critical 
temperature is investigated. Single-band models (BCS-like) with d-wave order parameter 
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display strong reduction of T c with increasing nonmagnetic impurity concentration 
Analogous effect on T c is produced by electron scattering by magnetic impurities in model 
with s-wave order parameter 0|TOj. In either case critical temperature vanish at certain im- 
purity concentration. However, the suppression of superconductivity induced by scattering 
from nonmagnetic impurities is inconsistent with existing experimental data [TT|, including 



the earlier ones, when samples were not pure enough. There are several different ways of 
possible explanation of weak reduction of critical temperature with increasing nonmagnetic 
impurity concentration in models with strongly anisotropic order parameter. The inclusion 
of anisotropy of scattering matrix element in the models with anisotropic order parameter 



12,13,14] can lead to weak impurity concentration dependence of T c . Relative arrangement 



of nodes of order parameter and those of impurity potential plays an important role here. It 



was demonstrated in |15| , that non-trivial topology of Fermi surface can strongly affect su- 
perconducting properties, and in particular inhibit reduction of T c with increasing impurity 
concentration. There is a separate class of models describing formation of order parameter 
anisotropy, namely, multi-band models. Although scattering by impurities was repeatedly 
investigated within the framework of these models [ |T^| , p75| , p7B| , [TT| , P~S| . [TD| , pD|] , a number of ques- 



tions concerning consideration of real crystal structure and symmetry properties of electron 
spectrum remains unsolved. 

The present article is concerned with effect of electron scattering by impurities on crit- 
ical temperature within the framework of previously proposed two-band model of an order 
parameter anisotropy ||. The width of one band is assumed to be much greater than that 
of another band. Hamiltonian of the model has the form 



Ho — ^2{e a ^ij + ttj)ai,a a j,t 
a iJ. a il a iT a «t 
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where the operators af^ a and create electrons on the i-th site with spin a in the wide 
and the narrow bands respectively; e a0 and e c0 are the energies of levels forming the wide 
and the narrow bands; and are the matrix elements of single-particle transitions 
between sites in the wide and the narrow bands; Wij is the matrix element of single-particle 
interband transitions. It is significant that since depend on the indices of different sites, 
in the k- representation the parameter W must be a function of k. It is assumed that 
superconductivity is caused by isotropic attraction in the narrow band U c < and there is 
isotropic effective interaction (either repulsive or attractive) between electrons in the wide 
band U a . 

The interaction of electrons with impurities is described by the following Hamiltonian 

Himp = Y, U aa (k, a; k', P)a^ a a k , ^ 

k,k',a,/3 

+ Y U cc (k,a;k',(3)c+ 

k,k',a,/3 

+ Y ( u ac(k, a; k',/3)a^ Q c k , )/3 

k,k',a,/3 

+h.c), (2) 

where 

i 

is the matrix element for electron scattering by impurities from the state (k',/3) in j band 
to the state (k, a) in i band = a, c), R/ is the position of l-th impurity . In general, 

uy(k, a; k', 0) = «.(k, k') + u£(k, k'))<W 

+^Jy(k,k')S i <x Q , / ? 

Here S; is the spin of magnetic impurity located at Hi, o a $ are the matrix elements of Pauli 
matrices, 

< M (k,k') = J ^t(r)u^\r)^\r)dh, 
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J^(k,k') = J4>(v)j(r)^ f \v)d 3 v, 

where u"( m )(r) is the potential of nonmagnetic (magnetic) impurity, J(r) is the exchange 
interaction of electron with magnetic impurity, ip^ is Bloch function of electron in i band. 

The impurity atoms will be assumed to be randomly distributed over the crystal lat- 
tice and the mean distance between them is much greater in comparison with interatomic 
distances. Let us consider short-range impurity potentials, so that the matrix elements 
u™j m \\s., k') and Jjj(k, k') do not depend on quasimomentum and are equal to u™^ 1 and 
Jij respectively. Thus, both intraband and interband electron scattering by impurities are 
isotropic. All calculations will be performed within Born approximation. 

Let us introduce the following Matsubara Green's functions: 

Gl j {r,r') = -i(T a i>a (r)af a (r')) , 

Gl j (T,T , ) = -i(Tc iia (T)cf a (T')), 

A, J (r,r / ) = -z(T a ha (r)c+ a (r')) , 

ioZpFgir, t') = -i (T aUr)al,{r')) , 

i< p F%{ T y) = -i(T < a (r) C +,(0), 

<pB±(T,T') = -i(T a+ a (r)cl p (r')) . 

The system of equations for these functions can be obtained by using equations of motion 
for the operators a itCe and Cj jQ ,. Then make transform to (k, ^-representation and perform 
Abrikosov-Gorkov procedure of averaging over impurity positions and, if impurities are mag- 
netic, over their spin orientations ||. As a result the system of equations for the narrow 
band Green's functions takes the following form: 

(iu c - e c (k)) G c (k, u) - W: a (k)D(k, u) 

+A c F+{k, w) + A ca B + (-k, -u) = 1, (3) 
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(iu a - e tt (k)) D(k, u) - W a * c (k)G c (k, u) 

+A a B+(-k, -u) + A ac F c +(k, u) = 0, (4) 

(iu* c + e c (k)) F c +(k, w) + ^ ca (k) J B+(-k, -w) 

+A c + G c (k,a;) + A+ J D(k,a;) = 0, (5) 

(iu* a + e a (k)) £? + (-k, -u) + ^ ac (k)F c + (k, 

+A a + J D(k^) + A+G c (k,c) = 0. (6) 

The system of equations for the wide band Green's functions can be obtained from (3-6) by 
interchanging of band indices a and c. 

Here new parameters are introduced, among which are renormalized frequencies u a , u c 
and anomalous means A a , A c appearing also in single-band models: 



A c 



= iuj — 


n (n) |<cl 2 + ^ (m) (KI 2 + 


\V CC \ 2 ) 


(2tt) 3 




X 


f d 3 kG c (k,uj) 






n {n) \<c\ 2 + n{m) (\ U ac\ 2 + 


Kc\ 2 ) 




(2vr) 3 




X 


J d 3 kG a (k,Lu), 




A c0 + 


n {n) \u n cc \ 2 + n^WuZf - 


M 2 ) 


(2tt) 3 




X 


J d 3 kF c {k,oo) 




+ 


n (n) \<c\ 2 + n{m H\ U ac\ 2 - 


Kd 2 ) 


(2tt) 3 




X 


/ d 3 kF a (k,u), 





(7) 



(8) 

where A c0 = —U c {c^cn), rS n "> and rS rn ^ are concentrations of nonmagnetic and magnetic 

impurities respectively |%| 2 = \j 2 j S{S-\- 1) = a,c), uj = (2n + 1)ttT. 

4 

In addition, electron scattering by impurities leads to renormalization of matrix element 
of hybridization W, as well as to appearance of interband anomalous self-energy A ca : 

(n) n n j. n (m)/..m..m i .. ? . \ 
W* = W* -I- u cc u ca ' n \ a cc a ca ' u cc u ca) 

(2tt) 3 
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, n^u n ca u n aa + n^juZuZ + v ca v aa ) 
(2vr) 3 

x J d 3 kG a (k,u;), (9) 

ca ~ (2^ 
x J d 3 kF c (k,uj) 

n (n)„n v n* i ^(rn) / m m* _ * \ 

I a ca a aa J[ ^ \ a ca a aa v ca v aa) 

(2vr) 3 

x | d 3 kF a (k,cu). (10) 

The expressions for cj a , A a , W ac and A ac can be obtained from (7-10) by interchanging 
of band indices a and c. 

The parameters A a0 and A c0 obey the self-consistency equations 

r d 3 k 

A a0 = -U a Tj2 J 7^r- 3 F a (k,u), (11) 



A 



c0 



r d 3 k 

WEf-^m^)- (12) 



d 3 k 

(2^ 

If impurities are absent in a superconductor (n^ = = 0), the system (4-7) is 

reduced to the Gorkov's system of equations for the functions G c and F c + with effective 
order parameter 

The essential feature of considered model is the strong anisotropy of the matrix element 
of single-particle interband hybridization W(k) which has lines of nodes along the diagonals 
of Brillouin zone for different types of symmetry of initial orbitals ||. As a result the 
effective order parameter (13) depends on quasimomentum in the case of nonzero interaction 
in the wide band. From Eqs. (11-13) it follows that if the interaction in the wide band is 
repulsive and its value is less than the critical one, at which supeconductivity in the system 
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is destroyed (U^ nt > U a > 0), then A a / A c < and the order parameter A c (k) changes its 

sign. Such a behaviour determines the main features of the model and makes it possible to 

explain some properties of HTSC compounds non-contradictory, e.g. strong anisotropy of 

the order parameter and experimentally observed s-like behaviour of — — characteristics of 

dv 

SIS' junctions 

The impurity concentration dependence of T c can be determined by solving the systems 
of equations (3-6) and (7-12) in the limit A a , A c — > at T — > T c . 

The chemical potential will be considered to be situated in the middle of the wide band 
with energy dispersion e a (k) = t a (cosk x + cosk y ). The location of the narrow band with 
energy dispersion e c (k) = e c o + t c (cosk x + cosk y ) relative to fi is defined by the parameter 
e C Q. The matrix element of hybridization, assuming different symmetries of initial orbitals, 
namely of s- and d- types, has the form jy(k) = Wo(cosk x — cosk y ). For numerical calcula- 
tions of impurity concentration dependencies of T c the following values of parameters were 
used: t a = 30, W = 5, U a = 27, U c = —8 (expressed in terms of t c ). 

The relative value of intraband and interband scattering matrix elements strongly de- 
pends on location of impurity atom in a unit cell. For example, these can be of the same 
order of magnitude, i.e. u aa oc u cc oc u ac . However, for substitutional impurities, scattering 
of electrons of one of the initial bands by impurities can overwhelm, while interband scat- 
tering turns out to be much weaker for symmetry reasons. It will be shown below that a 
behaviour of T c with increasing impurity concentration is different for various cases. 
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FIG. 1. The dependencies of the superconducting transition temperature T c on the parameter 
nu 2 corresponding to electron scattering by nonmagnetic impurities calculated for A) e c o = 0, B) 
e c0 = -3 and a) u = u c , u ac = u a = 0; b) u = u a , u ac = 0.5u a , u c = 0; c) u = u c , u ac = 0.5ii c , 
u a = 0; "s"- within the framework of single-band model with s-wave order parameter; "d"- within 
the framework of single-band model with <i-wave order parameter. 

In what follows the effect of electron scattering by nonmagnetic (Fig.l) and magnetic 
(Fig. 2) impurities on the critical temperature depending on the location of the center of 
the narrow band e C Q relative to the chemical potential fi is considered. The obtained T c (n) 
dependencies are compared with the results calculated within the framework of single-band 
models with s- and G?-wave order parameters. Here mean electron density of states of these 
models corresponds approximately to that of two-band model. 

Fig.l A gives the critical temperature T c as a function of parameter n^u 2 (where u = 
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u a a,Ucc,Uac), chemical potential being in the narrow band (e c o = 0). As can be seen in 
the figure, scattering in the narrow band has a strongest effect on the critical temperature. 
Including of interband scattering scarcely affects the result. This is supported by the fact 
that interband scattering together with intraband one in the wide band produces a weak 
effect on T c . Corresponding curve coincide with the one calulated within the framework of 
single-band model with s-wave order parameter very closely. 

If the center of the narrow band is far from the chemical potential (e c q = —3), the 
behaviour of T c with increasing impurity concentration is changed (Fig. IB). In this case 
interband scattering decreases the critical temperature while intraband scattering in the 
narrow band raises it in the region of small values of parameter rv^'u^,. The increase of 
T c at small impurity concentration can be explained by the fact that electron density of 
states increases in energy region near the chemical potential because of broadening of a 
peak originating from the van Hove singularity of the initial narrow band. The effect in this 
region of parameters overwhelms destruction of superconductivity caused by scattering by 
impurities. 

It is seen from Fig.l that though the order parameter of considered model changes its 
sign, nonmagnetic impurities suppress superconductivity much weaker than in single-band 
models with c?-wave order parameter. There is also no critical impurity concentration, at 
which the superconductivity in the system disappears. Intraband scattering in the wide 
band and interband scattering scarcely affect the critical temperature. 

A series of experiments aimed at investigating magnetic impurity concentration depen- 
dence of T c of HTSC compounds have been performed. As magnetic impurities Zn and Ni 
atoms substituting for copper sites in C11O2 plane were used P^j23|j2^j25[| . Experimental 



data indicates that these impurities influence critical temperature as impurity magnetic cen- 



ters rather than by changing carrier concentration p3[p6[ . It should be noted that while Zn 
substitution results in strong T c suppression in the copper oxides similar to that induced by 
magnetic impurities in conventional superconductors, Ni substitution was experimentally 
observed to affect T c weakly ||23|| . 
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FIG. 2. The dependencies of superconducting transition temperature T c on parameter nu 2 cor- 
responding to electron scattering by nonmagnetic (n) and magnetic (m) impurities calculated for 
A) e C Q = 0, B) e c o = —3 and a) u = -u™, ii" c = it" = 0; b) u = u™, u™ c = 0.5u™, u™ = 0; c) u = u™, 
< c = 0.5<, u n a = 0; a') « = <\ < c = < = 0; b') u = u™, u™ c = 0.5n™, u™ = 0; c') u = 
< = 0.5<, < = 0. 

Fig. 2 presents magnetic impurity concentration dependencies of T c . For comparison the 
T c (nu 2 ) curves for the case of nonmagnetic impurities are also shown in the figure. If the 
chemical potential is in the narrow band (Fig.2A), the quickest reduction of T c is caused by 
intraband scattering in it (a', c' curves). As in single-band models ||, the critical tempera- 
ture vanish at a certain impurity concentration. A comparison of b' and d curves shows that 
unusually weak effect on T c is produced by interband scattering from magnetic impurities. 
This fact as well as that intraband scattering scarcely affects the critical temperature can be 
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explained by small value of electron density of states in the energy region near the chemical 
potential and by the fact that superconductivity is due to attraction between electrons in the 
narrow band. The possibility of weak dependence of T c on magnetic impurity concentration 
by considering interband scattering was discussed earlier [^0] . 



If the center of the narrow band is far from the chemical potential (Fig.2B), then both 
intraband scattering in the narrow band and interband scattering results in the usual strong 
reduction of T c with increasing magnetic impurities concentration. 

The diagonalization of Hamiltonian of the considered model leads to the problem on 
anisotropic electron scattering by impurities in bands with highly anisotropic order param- 
eters which are roughly of s + d symmetry type. The symmetry of effective scattering 
matrix elements resulting from the diagonalization of the model Hamiltonian as well as that 
of the order parameter are determined by properties of electron spectrum. The degree of 
anisotropy of impurity potential depends on the symmetry of initial bands and on relative 
value of initial intraband and interband scattering matrix elements. The anisotropic compo- 
nent could also be introduced into initial impurity potentials, but this is beyond the scope 
of this article. 

Thus, within the framework of the two-band model of a superconductor with anisotropic 
order parameter based on symmetry properties of crystal lattice the impurity concentration 
dependence of the critical temperature was investigated. The performed calculations sug- 
gest that the results depend on both the location of the chemical potential relative to the 
narrow band and relative value of the intraband and interband impurity scattering matrix 
elements. It was found that the critical temperature is much less sensitive to the increase in 
concentration of nonmagnetic impurities than in one-band model with c?-wave order param- 



eter. As in the article |T7|], the conditions for T c increasing with impurity concentration were 
found. In accordance with standard theories strong dependence of the critical temperature 
on magnetic impurity concentration was found at some values of the model parameters. 

This work was supported by the RFBR Grant N02-02- 16925 and Landau Scholarship 
Grant. 
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